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ABSTRACT 

We have used a statistica l technique "Shuffle" (Bh avsar fc Line! 119881 : 
iBharadwai. Bhavsar fc Shethl 12004) in seven nearly two dimensional strips from the 
Sloan Digital Sky Survey Data Release Six (SDSS DR6) to test if the statistically 
significant length scale of filaments depends on luminosity, colour and morphology of 
galaxies. We find that although the average filamentarity depends on these galaxy 
properties, the statistically significant length scale of filaments does not depend on 
them. We compare it's measured values in SDSS against the predictions of ACDM 
N-body simulations and find that ACDM model is consistent with observations. The 
average filamentarity is known to be very sensitive to the bias parameter. Using ACDM 
N-body simulations we simulate mock galaxy distributions for SDSS NGP equatorial 
strip for different biases and test if the statistically significant length scale of filaments 
depends on bias. We find that statistically significant length scale of filaments is nearly 
independent of bias. This result is possibly related to the fact that statistically signif- 
icant length scale of filaments is nearly the same for different class of galaxies which 
are differently biased with respect to underlying dark matter distribution. The average 
filamentarity is also known to be dependent on the galaxy number density and size of 
the samples. We use ACDM dark matter N-body simulations to test if the statistically 
significant length scale of filaments depends on number density of galaxies and size 
of the samples. Our analysis shows that the statistically significant length scale of 
filaments very weakly depends on these factors. Finally we test the reliability of our 
method by applying it to controlled samples of segment Cox process and find that 
our method successfully recovers the length of the inputted segments. Summarizing 
these results we conclude that the statistically significant length scale of filaments is 
a robust measure of the galaxy distribution. 

Key words: methods: numerical - galaxies: statistics - cosmology: theory - cosmology: 
large scale structure of universe 



1 INTRODUCTION 

The fact that the galaxies appear to be distributed along 
filaments which are interconnected to form a web like struc- 
ture which is often reffered as the 'cosmic web' is one 
of the most striking visual feature in all the present and 
past re dshift surveys (e.g. , CfA , iGeller fc Huchra 1 19891: 
LCRS, IShectman etall Il99rj; 2d FGRS, IColles et al.l l200ll 
and SDSS. IStoughton et alj|2003 h 

The analysis of filamentary patterns in the galaxy dis- 
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tribution has a long history dating b ack to a few papers i n 
the late-seven tie s and mi d-eighties bv | joeveer et al.l (1 1978). 
lEinasto et alj (Il980l). IZel'dovich . Einasto fc Shandarin! 
dl982h . IShandarin fc Zeldovichl (|l983l ) and lEinasto et all 
(1984). Filaments are the most striking vis- 



(e.g. IGeller & Huchral Il989l. IShectman et al.l 


1996. 


Shandarin & Yess! 19981, Bharadwai et al.l 


2000. 


Miiller et al.l 2000, Basilakos, Plionis, & Rowan-Robinson 


2004 |Doroshkevich et all 


2004, 


Pimbblet, Drinkwater & Hawkrigd 


2004, 


Pimbblct & Drinkwater 2004, Pandcv & Bharadwai 2005). 



A review on a number of physically motivated and statistical 
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methods to define filaments is provi ded in IPimbbletl d2005l) 
The percolation a n alysis (eg. IShandarin fc Zeldovichl 
19831. lEinasto et al.l Il984l). th e genus statistics (eg 



Gott. Dickinson, fc Melott] 19861 ), the minimal spanning 



tree (e.g. iBarrow et al 



dlcke fc van de Weveaertl 




t he Voronoi tessellation 
Ivan de Weveaert fc Ickel 



19891). the Minkowski f unctio nal (eg. iMecke et al.l 1 1994 
Schmalzing fc Bucherj 1997 ) and th e 'Shapefinders' 



i Sahni. Sathvaprakash.fc Shandarinlll998h are some of the 
useful statistical tools introduced to qu antify the topology 
and geometry of the galaxy distribution. IStoica et al.l (|2005h 
propose to apply a marked point process to automatically 
deline ate filaments in the galaxy distribution. IColberg et ahl 
l|2005l ) studied the intercluster filamentary network in high 
resolution N-bod y simulations of structure formation in a 
ACDM Universe. fAragon-Calvo et al.l |2007h use Multiscale 
Morphology Filter technique to identify wall-like and 
filament-like s t ructu res in cosmological N-body simulations. 
ISousbie et al.l (|2008l ) propose a skeleton formalism to quan- 
tify t he filamentary struc ture in three dimensional density 
fields. [St pica et al. propose to apply an object point 

process t o objectively identify f ilame nts in galaxy redshift 
surveys. ISarkar fc Bharadwail (2009) propose the Local 
Dimension to locally quantify the shape of large scale 
structures in the neighbourhood of different galaxies in the 
Cosmic Web. 

The SDSS (|York et all |2000| ) is currently the 
largest galaxy redsh i ft su rvey. In an earlier work 
l|Pandev fc Bharadwail [2005) (hereafter Paper I) we 
have analysed the filamentarity in the equatorial strips of 
this survey. These strips are nearly two dimensional and 
we have projected the data onto a plane and analysed 
the resulting 2-D galaxy distribution. We find evidence 
for connectivity and filamentarity in excess of that of a 
random point distribution, indicating the existence of an 
interconnected network of filaments. We find that filaments 
are statistically sig nificant upto length scales 80fe~ 1 Mpc 
and not beyond (jPandev fc Bharadwail [2005). All the 
structures spanning length-scales larger than this length 
scale are the result of chance alignments. This is consistent 
with an earlier analysis by Bharadwai. Bhavsar fc Shethl 
(2004) where they show that in Las Campanas Redshift 
Survey (LCRS) the largest length-scale at which filaments 
are statistically significant is between 70 to 80ft _1 Mpc. 
Further we show that the average filamentarity of the 
galaxy distribution depends on various physical properties 
of galaxies such as luminosity, colour, morph o logy and 
star formation rate l|Pandev fc Bharadwail [200rj . 120081 ). It 
would be interesting to know if the statistically significant 
length scale of filaments also depends on different galaxy 
properties. In the present work we study if the statistically 
significant length scale of filaments depends on luminosity, 
colour and morphology of galaxies. 

Further it is also possible to measure the statistically 
significant length scale of filaments in mock galaxy distribu- 
tion extracted from N-body simulation. The ACDM model 
is currently believed to be the mini mal model which is con- 
sistent with most co s mological data |Efstathiou et al.ll20"oil; 
Percival etail l200ll; I Tegmark et all |2004 ISpergel et all 



and compare it with me asured values in the SDSS 
jPandev fc Bharadwail 120051 ). The N-body simulations pri- 
marily predict the clustering of the dark matter. This should 
be contrasted with the fact surveys reveal only the bright 
side of the matter distribution. Galaxy formation is a com- 
plicated process and the exact relation between the distri- 
bution of the galaxies and the dark matter is not well un- 
derstood. It is now generally accepted that the gala xies are 
a biased tracer of the da rk matter distribution (e.g.. iKaiserl 
ll984lMo fc Whitell 19961) and on large scales one expects the 
fluctuations in the galaxy and the dark matter distribution 
to be linearly related through the linear bias parameter b. 
Mock galaxy distributions with different bias can be sim- 
ulated following this assumption. The average filamentar- 
ity of the simulated galaxy dist ribution is found to be very 
sensitive to the bias parame ter l|Bharadwai fc Pandevll2004l ; 
iPandev fc Bharadwaill2007l ). It would be interesting to know 
how the statistically significant length scale of filaments de- 
pends on the bias parameter. In the present work we test if 
the statistically significant length scale of filaments depends 
on bias. 

Earlier we find that the average filamentarity is sen- 
sitive to_th£_area_aiid_g_alaxv_jrumber density of the sam- 
ples l|Pandev fc Bharadwail I2OO6I ) (hereafter Paper II) and 
this statistics can be used for a meaningful comparison be- 
tween two different galaxy samples only when they have the 
same volume (identical shape and size) and galaxy number 
density. In the present work we would also like to test if the 
statistically significant length scale of filaments depends on 
area and galaxy number density of the samples. 

Finally we simula te segment Cox process 
l|Pons-Borderia et al.lll999l ) with different segment lengths 
and test the efficiency of our method in measuring the 
statistically significant length scale of filaments in the mock 
samples drawn from these simulations. 

A brief outline of our paper follows. Section 2 describes 
the data and method of analysis, our results and conclusions 
are presented in Section 3. 



2 DATA AND METHOD OF ANALYSIS 



2.1 SDSS data 



20031 ; ISpergel et al.M2007l ; iKomatsu et al.| 120081 ). It would 
be interesting to measure the statistically significant length 
scale of filaments in N-body simulations of ACDM model 



The Sloan Digital Sky Survey (SDSS) l|York et al.ll2000D is 
a wide-field imaging and spectr oscopic survey of the sky 
using a dedicated 2.5 m telescope (|Gunn et al . 2006) with 3° 
field of view at Apache Point Observatory in southern New 
Mexico. The survey is carried out in five broad filters namely 
u, g, r, i a nd z covering the wavelength range from 30 00 to 
10000 A° (|Fukugita et al.lll996l ; Istoughton et al.ll2002D . 

Our present analysis is based on SDSS DR6 
|Adelman-McCarthv et alJlioolj ) galaxy redshift data. The 
SDSS DR6 includes 9583 deg 2 imaging and 7425 deg 2 of 
spectroscopy. The spectroscopic data includes 1,271,680 
spectra with 790, 860 galaxy redshift. Different selection al- 
gorithms are used for different categories of SDSS targets. 
We use the Main Galaxy Sample for the present work. The 
M ain Galaxy Sample target selection algorithm is detailed 
in IStrauss et al. | |2002| . The Main Galaxy Sample comprises 
of galaxies brighter than a limiting r band Petrosian magni- 
tude 17.77. We downloaded the data from Catalog Archive 
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Server (CAS) of SDSS DR6 by Structured Query Language 
(SQL) search. 

For the present analysis we use seven non-overlapping 
strips each spanning 90° in A and 2° in r\, lying entirely 
within the survey area of SDSS DR6. A and r/ a re sur- 
vey co-ordinates described in IStoughton et all <|2002h . These 
strips are identical in sky coverage as the ones used in 
IPandev fc Bharadwail |2006) and are shown in Figure 1 of 
that paper. For each strip we have extracted two different 
volume limited subsamples in binl and bin2 (Table 1 and 
Table 2). For each strip the number of galaxies in binl and 
bin2 are listed in Table 1. The absolute magnitude and red- 
shift limits for binl and bin2 are listed in Table 2. 

The analysis using "Shuffle" requires us to cut the entire 
survey area into squares and shuffle them around. The thick- 
ness of the resulting subsamples increases with redshift. For 
our analysis we have considered a smaller region of uniform 
thickness corresponding to the value at the lowest redshift. 
For each luminosity bin the area and number density are 
listed in Table 2. For each bin the galaxy number density 
varies slightly across the seven strips and the average den- 
sity along with the 1 — a variation is shown in Table 2. For 
all the subsamples, the thickness is much smaller than the 
other two dimensions and hence it is collapsed along the 
thickness resulting in a 2D distribution (Figure [!}. 

We study luminosity dependence of statistically signif- 
icant length scale of filaments by measuring and comparing 
it's value in binl and bin2. We label the galaxies in binl as 
Faint and galaxies in bin 2 as Bright. 

We also separately study the colour and morphology 
dependence in bin 2 (Table 1 and Table 2) by dividing 
the galaxies in this bin into Red/Blue and Early/Late type 
galaxies. The galaxies in bin 1 are not divided on the basis 
of their colour or morphology and hence for bin 1 no cut- 
offs in u — r colour and concentration index d are listed in 
Table 2. 

When testing for colour dependence, all the galaxies 
(Table 1) in the entire area of the bin 2 (Table 2) are classi- 
fied as either red or blue galaxies. The galax y u — r colours 
are k nown to have a bimodal distribution (jStrateva et al.l 
2001). In our analysis we determine a value (u — r) c for the 
colour such that it divides galaxies in bin 2 into equal num- 
ber of red (i.e. u — r > (u — r) c ) and blue (i.e. u — r < (u—r) c ) 
galaxies. As a consequence the number density of red and 
blue galaxies in bin 2 are exactly equal and have a value half 
that of the density shown in Table 2. 

The morphological classification was carried out us- 
ing the concentration index defined as Ci = rgo/rso where 
rgo and r$o are the radii containing 90% and 50% of 
the Petrosian flux respectively. This has been found to 
be o ne of the best paramet er to classify galaxy morphol- 
ogy (SShimasak u et all 1200 ll ). Ellipticals are expected to 
have a larger concentration indices than spirals. It was 
found that a ~ 3 .33 for a pure de-Vaucouleurs profile 
(Blan ton et al. 2001 ) while a ~ 2.3 for a pure exponential 
profile ( Strateva et al.ll200ll ) . In bin 2 we have chosen a cut- 
off Ci, c (Table 2) that partitions the galaxies into two equal 
halves, one predominantly ellipticals and the other spirals. 



2.2 N-body data 

We simulate the dark matter distribution using a Particle- 
Mesh (PM) N-body code. The simulations use 256 3 par- 
ticles on a 512 3 mesh, and they have a comoving volume 
^l.e/i^Mpc] 3 . We use (O ro0 ,«Ao,h) = (0.27,0.73,0.71) 
for the cosmological parameters along with a ACDM power 
spectrum w ith spectral index n B = 0.96 and normalization 
a 8 = 0.8 12 jKomatsu et al.ll2008T ). A "sharp cutoff" biasing, 
scheme |Cole et alj 1998?) was used to extract particles from 
the N-body simulations. Identifying these particles as galax- 
ies, we have galaxy distributions that are biased relative to 
the dark matter. The bias parameter b of each simulated 
galaxy sample was estimated using the ratio 



b = 



(i) 



where £ 9 (r) and £(r) are the galaxy and dark matter two- 
point correlation functions respectively. This ratio is found 
to be constant at length-scales r > 5/i _1 Mpc and we use the 
average value over 5 — 40/i _1 Mpc. We use this method to 
generate galaxy samples with bias values 1.2, 1.5 and 1.8. 
We also consider anti bias and galaxy distributions with 
bias b = 0.8 were generated by adding randomly distributed 
particles to the dark matter distribution. 

The peculiar velocity effects were included to produce 
galaxy distributions in redshift space. We have used three 
independent realisations of the N-body simulations, and for 
each value of bias we have extracted three different mock 
data-sets from each simulation. This gives us a total of nine 
simulated galaxy distribution for each bias values. The sim- 
ulated galaxy distributions (Figure QJ have the same area, 
thickness and number density as that of the uniform thick- 
ness volume limited subsample constructed out of equato- 
rial strip (145° < a < 236°, -1° < 8 < 1°) from Northern 
Galactic Cap (NGP) in Paper I. Each strip extends from 
235/i _1 Mpc to 571/i _1 Mpc comoving in the radial direc- 
tion, has uniform thickness of 8.2 ft -1 Mpc and contains 1936 
galaxies (Table 1 of Paper I). We choose to simulate this 
strip as we have already established in Paper I that filaments 
are statistically significant upto length scales 80 /i _1 Mpc and 
we want to compare this value with that measured from sim- 
ulations. The simulated data were analyzed in exactly the 
same way as the actual data. 

Earlier studies (Paper II) using N-body simulations 
show that the filamentarity in the galaxy distribution is sen- 
sitive to the area and galaxy number density of the samples. 
We want to test whether the statistically significant length 
scale of filaments depends on these factors. We extract re- 
gions of identical shape and size similar to the uniform thick- 
ness NGP strip in Paper I from dark matter ACDM N-body 
simulations. We label the mock SDSS NGP strip as density 
2 and area 3 in Table 3. Thus effectively the mock SDSS 
NGP strips, area 3 strips and density 2 strips (Table 3) rep- 
resent the same galaxy samples. We then prepare different 
sets of mock galaxy samples which have the same area as 
SDSS NGP strip but different number densities (density 1, 
density 3 in Table 3) and other which have the same number 
densities as SDSS NGP strip but cover different areas (area 
1 and area 2 in Table 3). For each sample described in Ta- 
ble 3 we extract three simulated strips from each simulation 
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Table 1. This shows the (A, if) range of the seven non-overlapping strips. For each strip the table shows the number of galaxies in volume 
limited subsamples bin 1 and bin 2 with absolute magnitude and redshift limits given in Table 2. 



Strip number 


A(°) 


v(°) 


bin 1 


bin 2 


1 


-50 < A < 40 


9 < ri < 11 


845 


1519 


2 


-50 < A < 40 


11 < V < 13 


743 


1253 


3 


-60 < A < 30 


13 < T) < 15 


624 


1153 


4 


-60 < A < 30 


15 < r) < 17 


718 


1188 


5 


-50 < A < 40 


21.5 < r? < 23.5 


774 


1381 


6 


-50 < A < 40 


24 < r? < 26 


799 


1338 


7 


-50 < A < 40 


26 < r] < 28 


781 


1239 



Table 2. This shows the absolute magnitude and redshift limits for the different volume limited subsamples analyzed. The area and the 
average galaxy number density with 1 — a variations from the 7 strips are also shown. The last two column shows the value of the galaxy 
colour (u — r) c and the value of the concentration index Ci jC used to divide the data into equal numbers of red/blue and elliptical/spiral 
galaxies respectively. Galaxies in bin 1 are not divided and only galaxies in bin2 are divided to obtain different samples for red/blue and 
elliptical/spiral galaxies. 



bin 


Absolute Magnitude range 


Redshift range 


Area [ 10 4 /i" 2 Mpc 2 ] 


Density [10 -2 /i 2 Mpc~ 2 ] 


(u - r) c c ijC 


bin 1 


-19 > M r > -20 


0.028014 < z < 0.075635 


3.35 


2.25 ± 0.19 




bin 2 


-20 > M r > -21 


0.043657 < z < 0.114635 


7.47 


1.73 ± 0.15 


2.38 2.7 
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Figure 1. This figure exhibits how Shuffle works. The left 
panel shows the galaxy distribution in a simulated NGP strip 
from dark matter ACDM N-body simulations and the right panel 
shows a shuffled realization generated from the same data using 
L = 90/i _1 Mpc. A 90/i _1 Mpc X 90/i _1 Mpc grid was placed on 
the mock NGP data, and square blocks lying fully inside the sur- 
vey region were randomly shuffled around to generate the data 
shown in the right panel. This process destroys all coherent struc- 
tures spanning length-scales larger than L = 90h~ x Mpc in the 
actual data and filaments larger than this in the shuffled data 
arise from chance alignments. Among the filaments seen in the 
right panel, those which run across the block boundaries have 
formed purely from chance alignments. 

giving us nine simulated slices. All these simulated slices are 
analyzed in exactly the same way as the actual data. 

2.3 Segment Cox Process 

We use the segment Cox process (|Pons-Borderia et al.lll999l ) 
to test the reliability of our method in determining the sta- 
tistically significant length scale of filaments. This is a con- 
trolled point process in which the segments of length I are 
scattered with random positions and orientations within a 
cube. Points are then randomly distributed on these seg- 
ments. The length density of the system of segments is 
L v — X 3 l, where X s is the mean number of segments per 
unit volume. If A; is the mean number of points on a seg- 
ment per unit length, then the intensity A of the resulting 
point process is, A = A; L v = A; \ s I. 




° 0.2 0.4 0.6 0.2 0.4 0.6 0.8 



Filling Factor (FF) 

Figure 2. This shows the Average Filamentarity as a function 
of Filling Factor (FF) for SDSS galaxies in strip 1 (Table 1) to- 
gether with the results for the shuffled data for two values of 
L shown in the figure. Results for different class of galaxies are 
shown in different panels and indicated in each panel. In all the 
panels shuffling with L = 10h~ 1 Mpc causes a large drop in the 
Average Filamentarity showing the statistical significance of the 
filamentarity at this length-scale. The data is within the 1 — a 
error bars of the shuffled realizations for L = 90h~ 1 Mpc (top left 
panel), L = 70/i _1 Mpc (top right panel) and L = 80/i _1 Mpc 
(bottom left and right panel). The filamentarity is statistically 
significant up to Lmax = 80,60 and 70/i -1 Mpc in these cases 
where the actual data lies just above the the 1 — a error bars and 
are not shown in this figure. For all the larger values of L, the 
data remains within the 1 — a error bars of the shuffled realiza- 
tions indicating that the filaments are not statistically significant 
beyond L MAX . 

For each of the two different segment length (I — 
20/i~ 1 Mpc and I = 80 fe _1 Mpc) we simulate three realiza- 
tions of segment Cox process inside a cube of side length 
L — 921.6 h~ Mpc. We choose the same box size as the N- 



Statistically significant length scale of filaments as a robust measure of galaxy distribution 5 



Table 3. This shows the redshift limits, area and density for the different mock subsamples analyzed to test the effect of area and density 
on the statistically significant length scale of filaments. 



sample Redshift range Area [ 10 4 h~ 2 Mpc 2 ] Density [10~ 2 h 2 Mpc -2 ] 

area 1 0.01 < z < 0.1 6.66 9.09 

area 2 0.05 < z < 0.15 13.06 9.09 

area 3 (NGP strip) 0.08 < z < 0.2 21.3 9.09 



density 1 0.08 < z < 0.2 21.3 4.54 

density 2 (NGP strip) 0.08 < z < 0.2 21.3 9.09 

density 3 0.08 < z < 0.2 21.3 18.17 



O 



X o 





Bright- 
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Red _ 
******** 


Blue _ 




Early _ 


Late _ 



30 60 90 120 30 60 90 120 150 



L (Mpc/h) 

Figure 3. This shows x 2 / u a * different shuffling lengths for 
different class of SDSS galaxies as indicated in each panel. The 
black dotted line indicates x 2 / u = 1- 
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Filling Factor (FF) 

Figure 4. This shows the results similar to Figure 2 but for a 
simulated SDSS NGP strip with different bias values as indicated 
in each panel. 



process for two different segment lengths are listed in Ta- 
ble 4. We keep the intensity of the resulting point process 
same in two cases by setting A; fixed and adjusting the value 
of A 3 with the value of I. Thus we have simulated the segment 
Cox process for two different values of / and for each val- 
ues of I the point process is simulated inside three different 
boxes. We extract three different simulated strips from each 
of the three boxes giving us a total nine simulated galaxy 
distribution for each I values. The simulated strips have the 
same area, thickness and the number density as that of the 
uniform thickness SDSS NGP strip used in Paper I. The 
simulated strips are analyzed in exactly the same way as 
the actual data. 

The simulated 2D strips are constructed out of 3D sim- 
ulations of segment Cox process by slicing and projecting 
and this fact destroys many segments and shortens many 
others. We aim to quantify this effect by generating pla- 
nar (2D) segment Cox process where all the segments lie on 
the same plane. We simulate the 2D segment Cox process 
over a planar region which has area and geometry identical 
to the uniform thickness SDSS NGP strip used in Paper I. 
The segments of length I are scattered with random posi- 
tions and orientations within this planar region. Points are 
then randomly distributed on these segments. The intensity 
A of the resulting point process is, A = A; A s Z where A s is 
the mean number of segments per unit area, A; is the mean 
number of points on a segment per unit length and I is the 
length of the segments. We simulate 9 such realizations of 
2D segment Cox process for each of the two segment lengths 
(/ = 20/i" 1 Mpc and I = 80/i _1 Mpc). The different param- 
eters used to simulate the 2D segment Cox process for two 
different segment lengths are listed in Table 5. We now ran- 
domly extract exactly same number of points as there are in 
the uniform thickness SDSS NGP strip used in Paper I from 
each of the 9 realizations of 2D segment Cox process. This 
gives us a total nine simulated galaxy distribution for each 
segment lengths (Table 5). The simulated strips have same 
area, geometry and number density as the SDSS NGP strip. 
We analyze these simulated strips exactly the same way as 
the actual data. 



body simulations used in this paper to enable us to extract 
regions of identical shape and size similar to the uniform 
thickness NGP strip used in Paper I. Two different segment 
lengths are chosen to represent shorter and longer filaments. 
The different parameters used to simulate the segment Cox 



2.4 Method of Analysis 

All the strips that we have analyzed are nearly two di- 
mensional. The strips were all collapsed along the thickness 
(the smallest dimension) to produce 2D gal axy distributions. 
We use the 2D "Shapefinder" statistic l|Bharadwai et all 
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Table 4. This shows different parameters used to simulate 3D segment Cox process. 



1 


A a 


A; 


A 


(Segment length) 


(Mean no. of segments per unit volume) 


(Mean no. of points on a segment per unit length) 


(Intensity) 


[h- 1 Mpc] 


[h 3 Mpc~ 3 ] 


[fcMpr 1 ] 


[h 3 Mpc" 3 ] 


20 


0.001278 


0.8 


0.020441 


80 


0.000319 


0.8 


0.020441 



Table 5. This shows different parameters used to simulate 2D segment Cox process. 



1 


A s 


A; 


A 


(Segment length) 


(Mean no. of segments per unit area) 


(Mean no. of points on a segment per unit length) 


(Intensity) 


[ft- 1 Mpc] 


[h 2 Mpc' 2 ] 


[hMpc- 1 ] 


[h 2 Mpc" 2 ] 


20 


0.037612 


0.8 


0.601791 


80 


0.009403 


0.8 


0.601791 



2000) to quantify the average filamentarity of the pat- 
terns in the resulting galaxy distribution. A detailed 
discussion is presented in Paper II , and we present 
only the salient features here. The read er is referred to 
ISahni. Sathvaprakash.fe Shandarinl (|l998h for a discussion 
of Shapefinders in three dimensions. 

The galaxy distribution is represented as a set of Is 
on a 2-D rectangular grid of spacing 1 h~ Mpc x 1 ft -1 Mpc, 
empty cells are assigned a value and filled cells are as- 
signed a value 1. We identify connected filled cells using 
the 'Friends-of-Friend' (FOF) algorithm. The filamentarity 
of each cluster is quantified using the Shapefinder T defined 



T = 



(P 2 - 16g) 
(P - Al) 2 



(2) 



where P and S are respectively the perimeter and the area 
of the cluster, and I is the grid spacing. The Shapefinder T 
has values and 1 for a square and filament respectively, 
and it assumes intermediate values as a square is deformed 
to a filament. We use the average filamentarity 



F 2 



E 



Sfj- i 



(3) 



to asses the overall filamentarity of the clusters in the galaxy 
distribution. 

The distribution of Is corresponding to the galaxies is 
sparse. Only ~ 1% of the cells contain galaxies and there 
are very few filled cells which are interconnected. As a con- 
sequence FOF fails to identify the large coherent struc- 
tures which corresponds to filaments in the galaxy distribu- 
tion. We overcome this by successively coarse-graining the 
galaxy distribution. In each iteration of coarse- graining all 
the empty cells adjacent to a filled cell (i.e. cells at the 4 
sides and 4 corners of a filled cell) are assigned a value 1 
. This causes clusters to grow, first because of the growth 
of individual filled cells, and then by the merger of adja- 
cent clusters as they overlap. Coherent structures extending 
across progressively larger length-scales are identified in con- 
secutive iterations of coarse-graining. So as not to restrict 
our analysis to an arbitrarily chosen level of coarse-graining, 
we study the average filamentarity after each iteration of 
coarse-graining. The filling factor FF quantifies the fraction 



of cells that are filled and its value increases from ~ 0.01 and 
approaches 1 as the coarse-graining proceeds. We study the 
average filamentarity F2 as a function of the filling factor 
FF (Figure |2| as a quantitative measure of the filamentar- 
ity at different levels of coarse-graining. The values of FF 
corresponding to a particular level of coarse-graining shows 
a slight variation from strip to strip. In order to combine 
and compare the results from different strips, for each strip 
we have interpolated F% to 7 values of FF at an uniform 
spacing of 0.1 over the interval 0.05 to 0.65. Coarse-graining 
beyond FF ~ 0.65 washes away the filaments and hence we 
do not include this range for our analysis. 

We use a statistical technique "Shuffle" to determine 
the largest length-scale at which the filamentarity is statis- 
tically significant. A grid with squares blocks of side L is 
superposed on the original data slice (Figure [1]). Blocks of 
data which lie entirely within the slice are then randomly 
interchanged, with rotation, repeatedly, to form a new shuf- 
fled slice. The shuffling process eliminates coherent features 
in the original data on scales larger than L, keeping cluster- 
ing at scales below L nearly identical to the original data. All 
the structures spanning length-scales greater than L that ex- 
ist in the shuffled slices are the result of chance alignments. 
At a fixed value of L, the average filamentarity in the orig- 
inal sample will be larger than in the shuffled data only if 
the actual data has more filaments spanning length-scales 
larger than L, than that expected from chance alignments. 
The largest value of L, Lmax, for which the average filamen- 
tarity of the shuffled slices is less than the average filamen- 
tarity of the actual data gives us the largest length-scale at 
which the filamentarity is statistically significant. Filaments 
spanning length-scales larger than Lmax arise purely from 
chance alignments. 

For each value of L we generated 24 different realiza- 
tion of the shuffled slices. To ensure that the edges of the 
blocks which are shuffled around do not cut the actual fil- 
amentary pattern at exactly the same place in all the real- 
izations of the shuffled data, we randomly shifted the origin 
of the grid used to define the blocks. The values of FF and 
F2 in the 24 realizations differ from one another and from 
the actual data at the same stage of coarse-graining. So as 
to be able to quantitatively compare the shuffled realiza- 
tions with the actual data, we interpolate the values of F2 
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in the shuffled realization at the same values of FF as inter- 
polated for the actual data. The mean ^[Shuffled] and the 
variance (AF 2 [Shuffled]) 2 of the average filamentarity was 
determined for the shuffled data at each value of FF using 
the 24 realizations. The difference between the filamentarity 
of the shuffled data and the actual data was quantified using 
the reduced x 2 P er degree of freedom 

= J_ V (^[Actual] [Shuffled])) 2 
v N p 2^ (AF 2 [Shuffled])2 ( > 

a— 1 

where the sum is over different values of the filling factor 
FF. 



3 RESULTS AND CONCLUSIONS 

We use "Shuffle" to determine the statistically significant 
length scale of filaments in the seven SDSS strips in mag- 
nitude bin 1 and bin 2 (Table 1 and Table 2). We label 
bin 1 as Faint and bin 2 as Bright. The top two panels of 
Figure [2] shows the results of "Shuffle" on strip 1 from bin 
1 and bin 2 respectively. In top two panels of Figure we 
see that shuffling the data with L = 10 causes a large drop 
in the average filamentarity (Fb). The average filamentar- 
ity increases as the value of L increases and the value of 
average filamentarity of the unshuffled slice comes within 
1 — a error bars when the data is shuffled with L = 90 and 
L — 70 for Bright and Faint galaxies respectively. Difference 
between the average filamentarity of shuffled and unshuffled 
slice is quantified by \ 2 l v an d the top two panels of Figure 
shows x 2 l v as a function of L. In top two panels of Fig- 
ure |3] we see that x 2 l v approaches 1 at L = 90 for Bright 
and L = 70 for Faint galaxies. So in strip 1 from binl and 
bin 2 the filaments are statistically significant upto length 
scales 80/i -1 Mpc and 60fo -1 Mpc respectively. The results 
from the other strips are not shown here. Averaging the re- 
sults from all the nine strips we find that the filaments are 
statistically significant upto length scales 84 ± 16/i _1 Mpc 
and 71 ± 16 h~ 1 Mpc in bin 1 and bin 2 respectively. 

The right and left middle panels of Figure [2] shows the 
results of "Shuffle" for Red and Blue galaxies (Table 2) from 
strip 1 in bin 2. The right and left middle panels of Figure 
[2] shows that a large drop in the average filamentarity is ob- 
served when the data is shuffled with L = 10. In both the 
panels we see that shuffling the data with L — 80 increase 
the average filamentarity and brings it back in agreement 
with the actual data. In two middle panels of Figure|3]we see 
that in both class of galaxies X 2 l v approaches 1 at L = 80 es- 
tablishing that for both Red and Blue galaxies the filaments 
are statistically significant upto length scales 70/i -1 Mpc. 
Here we have not shown the results from the other strips. 
Combining results from all the nine strips we find that 
the filaments are statistically significant upto length scales 
77±10 /i _1 Mpc and 71±10 /i _1 Mpc for Red and Blue galax- 
ies respectively. 

We show the results of "Shuffle" for the Early and Late 
type galaxies (Table 2) in the right and left bottom panels 
of Figure [2] Here also we see that the average filamentar- 
ity decreases when the data is shuffled with L — 10 and it 
increases and matches with the actual data when the data 
is shuffled with L = 80. The x 2 l v approaches 1 at L = 80 



for both Early and Late type galaxies as shown in two bot- 
tom panels of Figure [3] This establishes that for both Early 
and Late type galaxies the filaments are statistically signifi- 
cant upto length scales 70 h~ 1 Mpc and not beyond. We have 
shown here only the results for strip 1 in bin 2. Combining 
the results from all the strips the filaments are found to be 
statistically significant upto length scales 75 ± 10/i _1 Mpc 
and 74 ± 17/i -1 Mpc for Early and Late type galaxies re- 
spectively. 

It is to be noted in middle and bottom two panels of 
Figure [5] that the Red and Early type galaxies have higher 
degree of filamentarity as compared to Blue and Late type 
galaxies at smaller values of fill i ng fac tor. We reported this 
effect in iPandev fc Bharadwail £2006). It is interesting to 
note that although the average filamentarity (7*2 ) depends 
on luminosity, colour and morphology of galaxies, the statis- 
tically significant length scale of filaments does not depend 
on these galaxy properties. 

In Figure|3]we show the results of "Shuffle" in simulated 
SDSS NGP strips with different biases. The top left panel 
show the result for the dark matter distribution from ACDM 
N-body simulations. In this case galaxies are assumed to ex- 
actly trace the dark matter and the bias parameter 6=1. 
We see that shuffling the simulated data with L = 10 causes 
a large drop in the average filamentarity of the simulated 
galaxy distribution. With increasing L the value of average 
filamentarity of the shuffled slice slowly approaches the val- 
ues corresponding to the unshuffled data. It is found that at 
L — 90 the two are within 1 — a error bars. The top right 
panel shows the result for b=1.2. We see a very similar re- 
sult and the unshuffled and shuffled data agrees well when 
the data is shuffled with L = 90. In top left and right panels 
of Figure [5] we see that X 2 l v approaches 1 at L — 90 for 
both 6=1 and 6 = 1.2 and hence filaments are statistically 
significant upto length scales 80 /i -1 Mpc for both cases. The 
results are shown for one single mock SDSS strip. Averaging 
the results from all the nine strips it is found that for b=l 
and b=1.2 filaments are statistically significant upto length 
scales 93±16ft~ 1 Mpc and 92 ± 15 h~ 1 Mpc respectively. The 
bottom left panel of Figure 2] shows the result of shuffling on 
simulated SDSS strips for a bias b=0.8. As usual shuffling 
the data with L = 10 causes a drop in average filamentar- 
ity of the shuffled slice but the drop is smaller as compared 
to other bias values. As L increases the average filamentar- 
ity slowly grows and approaches the values corresponding to 
the unshuffled data at L = 90. In the bottom right panel of 
Figure [4] results are shown for a high bias b=1.5. Shuffling 
shows similar results for b=1.5 but a relatively large drop is 
seen when the data is shuffled with L = 10. With increasing 
L the average filamentarity grows relatively faster as com- 
pared to b=0.8 and finally the average filamentarity in the 
shuffled data levels up with the unshuffled data when L = 80 
is used for shuffling. In the bottom left and right panels of 
Figure [5] it is found that x 2 l v approaches 1 at L = 90 and 
L — 80 indicating that filaments are statistically significant 
upto length scales 80 h~ 1 Mpc and 70 /i _1 Mpc for b=0.8 and 
b=1.5 respectively. These are the results shown for a single 
mock SDSS strip and by combining the results from the 
other strips we get that the filaments are statistically signif- 
icant upto length scales 91 ± 13 /i _1 Mpc and 90± 15 /i -1 Mpc 
for 6 = 0.8 and 6 = 1.5 respectively. 

It is to be noted in Figure [4] that as the bias 6 is in- 
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Figure 5. This shows \ 2 Jv at different shuffling lengths for a 
simulated SDSS NGP strip with different bias values as indicated 
in each panel. The black dotted line indicates x 2 / u = 1- 
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Figure 6. This shows the statistically significant length scale of 
filaments as a function of bias b. The error bars are 1 — a error 
bars measured from nine mock galaxy catalogues. 

creased the average filamenta rity (F%) increases at smalle r 
FF. We reported this effect in lBharadwai fc Pandevl l|2004) . 
But interestingly the shuffling length at which the avergae 
filamentarity of the unshuflled and shuffled data agrees does 
not depend on bias. Figure [6] shows the statistically signif- 
icant length scale of filaments as a function of bias. This 
Figure shows that the statistically significant length scale 
of filaments is nearly independent of bias. The mean value 
is 90 h~ Mpc for the bias values b = 0.8 — 1.5. It tends to 
decrease for a high value of bias b — 1.8. This could be pos- 
sibly related to the fact that at very high values of bias we 
preferentially identify only very high density regions which 
are mostly related to cluster like structures rather than fil- 
aments. 

Earlier in Paper I we established that the filaments 
are statistically significant upto length scales of 80 ft -1 Mpc. 
This measured value lies well within 1 — a error bars of the 
measured values for mock galaxy distribution from ACDM 
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Figure 7. This shows the statistically significant length scale of 
filaments as a function area of the samples described in Table 3. 
The error bars are 1 — a error bars measured from nine mock 
galaxy catalogues. 

o 

CN I 1 1 1 1 



O 




o 

CO 



o 

1^ 1 

o I I I I 

5x10~ 3 0.01 0.015 0.02 

Density ((Mpc/h)" 2 ) 

Figure 8. This shows the statistically significant length scale 
of filaments as a function galaxy number density of the samples 
described in Table 3. The error bars are 1 — <r error bars measured 
from nine mock galaxy catalogues. 

N-body simulations. We conclude that the ACDM is consis- 
tent with SDSS observations. 

The galaxy samples constructed over different magni- 
tude ranges have different area and galaxy number density 
(Table 2). Earlier analysis in Paper II shows that the aver- 
age filamentarity of the galaxy distribution depends on the 
area and galaxy number density of the galaxy samples and 
filamentarity between different galaxy samples can only be 
compared when the samples have same area, geometry and 
number density of galaxies. So it is important to know if 
the statistically significant length scale of filaments depends 
on these factors before we can compare it's measured val- 
ues among various galaxy samples having different areas and 
galaxy number densities. 
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We prepare two different set of mock galaxy samples 
namely area 1 and area 2 (Table 3) which have the same 
galaxy number density as the SDSS NGP strip (area 3) but 
cover different redshift ranges and hence have different areas 
as listed in Table 3. We extract three mock galaxy samples 
from each of the ACDM dark matter simulation for each set 
of galaxy samples area 1 and area 2. We determine the sta- 
tistically significant length scale of filaments for these galaxy 
samples. The analysis for the mock SDSS NGP strip which 
we name as area 3 in Table 3 are already done. We plot the 
statistically significant length scale of filaments as a func- 
tion of the area of the galaxy samples in Figure [7] From this 
figure we see that the samples of smaller area tend to have 
a smaller length scale of filaments. But this trend is weak 
and the size of the error bars are quite large. We conclude 
that the statistically significant length scale of filaments are 
nearly independent of the area of the samples. 

We next prepare two different set of mock galaxy sam- 
ples namely density 1 and density 3 (Table 3) which have 
the same area as the SDSS NGP strip (density 2) but have 
different number densities listed in Table 3. We extract three 
mock galaxy samples from each of the ACDM dark matter 
simulation for each set of galaxy samples density 1 and den- 
sity 3. We determine the statistically significant length scale 
of filaments for these galaxy samples. The analysis for the 
mock SDSS NGP strip which we name as density 2 in Ta- 
ble 3 are already done. We plot the statistically significant 
length scale of filaments as a function of the number den- 
sity of galaxies of the samples in Figure [8] We see a trend 
of low density samples to have a smaller length scale of fil- 
aments. But the error bars are quite large and the depen- 
dence is weak. We conclude that the statistically significant 
length scale of filaments weakly depends on the number den- 
sity of galaxies. In a recent work with SDSS LRG samples 
|Pandev et alj2009f ) having much larger area and lower den- 
sity we find that the statistically significant length scale of 
filaments does not change much as compared to it's value 
measured in SDSS Main galaxy samples. Thus the statisti- 
cally significant length scale of filaments is a more robust 
statistics than the average filamentarity. It does not depend 
on properties (size and number density) of galaxy samples 
, and physical properties of galaxies. It emerges as a robust 
measure of the large scale structures. The statistically sig- 
nificant length scale of filaments could be also thought of as 
an indicator of the scale beyond which the galaxy distribu- 
tion become homogeneous because it tells us the length scale 
beyond which no coherent structures exist in th e galaxy dis- 
tribu tion. Multifractal analysis in SPSS DR1 rtYaday et al.l 
120051 ) and recently in SDSS DR6 (|Sarkar et aUl2009h show 
that the galaxy distribution becomes homogeneou s at a 
length scale between 60 to 70/i _1 Mpc. lYadav et aD (|2005l ) 
considered the transition to homogeneity in ACDM simula- 
tions with different bias 6=1, 1.6, 2 and find that irrespec- 
tive of the bias values the distribution become homogeneous 
at a length scale between 60 to 70ft" 1 Mpc. These results 
are consistent with our findings. 

Finally we test the reliability of our method by apply- 
ing it to different controlled samples (Table 4) of 3D seg- 
ment Cox process. We use segments of two different length 
I = 20/i _1 Mpc and I = 80ft _1 Mpc for simulating the 3D 
segment Cox process. Two different sets of 2D mock sam- 
ples drawn from these 3D simulations (Table 4) are sepa- 



rately analyzed. The results are shown in top two panels of 
Figure [10] and Figure 1111 We see in top left panel of Fig- 
ure [TO] that filamentarity decrease by a very small amount 
when the mock samples having I = 20 ft -1 Mpc are shuf- 
fled with L — 10 and there are virtually no drop in fila- 
mentarity when the data is shuffled with L = 20. The top 
right panel of the same figure shows the result of "Shuffle" 
on cotrolled samples of segment Cox process with segment 
length I = 80 ft" 1 Mpc. We see a significant drop in the av- 
erage filamentarity when the data is shuffled with L = 10 
and the average filamentarity finally saturates to the val- 
ues corresponding to the unshuffled data when L = 60 is 
used for shuffling. The \ 2 l v as a function of L for these 
two cases are shown in top two panels of Figure 1111 We see 
that x 2 / u approaches 1 at L = 20 and L — 60 in these two 
cases respectively. This indicates that for these 2D mock 
samples having I = 20/i~ 1 Mpc and I = 80 h~ Mpc, the fil- 
aments are found to be statistically significant upto length 
scales 10/i _1 Mpc and 50 ft -1 Mpc respectively. Here the re- 
sults are shown only for a single simulated mock sample. 
By averaging the results from all the nine simulated mock 
samples we find that in the above two cases the filaments 
are statistically significant upto length scales 18±7/i~ 1 Mpc 
and 52 ± 11 ft - 1 Mpc respectively. We note that the values 
of statistically significant length scale of filaments recovered 
by our method for straight segments simulated here are a 
little shorter but close to the original values of I used in the 
simulations. The simulations are carried out in 3D boxes 
and our mock samples are the 2D projections of nearly two 
dimensional samples drawn from these boxes. This possibly 
destroys and shortens a lot of segments and the effects are 
expected to be more prominent for longer segments as re- 
flected in our results. This effect indicates that also for the 
2D real galaxy samples from SDSS which are constructed 
out of a 3D galaxy distribution, Shuffle possibly underesti- 
mates the values of statistically significant length scale of 
the real 3D filaments. 

The effects of slicing and projections can be avoided if 
we simulate the segment Cox process in 2D as in this case 
the segments lie on the same plane (Figure [§} . We apply 
Shuffle to a set of mock 2D strips constructed out of sim- 
ulations of 2D segment Cox process (Table 5). We use two 
different segment lengths I — 20 /i _1 Mpc and I — 80 /t -1 Mpc 
to simulate the 2D segment Cox process. We use the same 
segment lengths as used in the simulations of 3D segment 
Cox process so as to test the effect of slicing and projection. 
We show the results in bottom two panels of Figure [10] and 
Figure 1111 Bottom left panel of Figure [TO] shows that when 
the mock samples with / = 20/i _1 Mpc are shuffled with 
L — 10 the average filamentarity of the shuffled data falls 
below that of the unshuffled data. The average filamentarity 
of the shuffled data saturates to the unshuffled data when 
the data is shuffled with L — 40. The bottom right panel of 
Figure [TO] shows a significant drop in the average filamentar- 
ity when the mock samples with I = 80 h~ Mpc are shuffled 
with L — 10. The average filamentarity of the shuffled data 
increases slowly with increasing shuffling length and finally 
saturates to the unshuffled data at shuffling length L = 90. 
We see in bottom two panels of Figure[TTJthe x 2 /v as a func- 
tion of L approaches 1 at L = 40 and L — 90 in these two 
cases respectively. This indicates that for these controlled 
samples of 2D segment Cox process having I = 20 ft -1 Mpc 
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Figure 9. This shows a simulation of 2D segment Cox pro- 
cess with two different segment length I = 20 h~ x Mpc and 
£ = 80fe _1 Mpc (as indicated in each panel) over a region which 
has identical geometry to the SDSS NGP equatorial strip. These 
simulated strips have same values of I and A; listed in Table 5 but 
different values of A s are chosen so as to lower the density of the 
segments and make them visibly distinguishable in this plot. 



and / = 80 h~ Mpc, the filaments are statistically signifi- 
cant upto length scales 30/i -1 Mpc and 80/i -1 Mpc respec- 
tively. We have shown the results for a single simulated 
strip in these plots. We combine the results from all the 
nine simulated strips and find that the filaments are sta- 
tistically significant upto length scales 25 ± 8/i _1 Mpc and 
82 ± 24/i~ 1 Mpc for simulations of 2D segment Cox process 
with segment length I = 20h~ 1 Mpc and / = 80/i _1 Mpc 
respectively. We note that our method can successfully re- 
cover better the length of inputted segments when the mock 
samples are drawn from simulations of 2D segment Cox pro- 
cess where the effects of slicing and projections are absent. 
There would be another effect due to the boundaries of the 
samples which again shortens the length of some segments 
near the boundaries but this effect is not so important in our 
analysis as the extent of the samples are much larger than 
the length of the inputted segments. This analysis suggets 
that the actual length scale of the real 3D filaments can be 
recovered with Shuffle if the analysis is extended in 3D. We 
propose to take up this in a future work. 

In conclusion we note that although the average fila- 
mentarity of the galaxy distribution depends on different 
galaxy properties, the statistically significant length scale of 
filaments does not depend on luminosity, colour and mor- 
phology of galaxies. We find that the measured values of 
the statistically significant length scale of filaments in SDSS 
are consistent with that measured from ACDM simulations. 
We considered the ACDM model with different values of 
bias and find that the statistically significant length scale 
of filaments is nearly independent of bias. Different class 
of galaxies are differently biased with respect to underly- 
ing dark matter distribution and this result is possibly re- 
lated to the fact that statistically significant length scale of 
filaments is nearly the same for different class of galaxies. 
Unlike the average filamentarity, the statistically significant 
length scale of filaments is nearly independent of the area 
and galaxy number density of the samples. This establishes 
the statistically significant length scale of filaments as a ro- 
bust statistics of galaxy distribution. 
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Figure 10. This shows the Average Filamentarity as a function 
of Filling Factor (FF) for a mock SDSS NGP strip constructed 
out of 3D simulation of segment Cox process (Table 4) and 2D 
simulation of segment Cox process (Table 5) together with the 
results for the shuffled data for two values of L shown in the 
figure. The 2D simulations of segment Cox process contain all 
the segments in the same plane and there are no effects because 
of the slicing process or the projections which arises when samples 
are constructed out of 3D simulations of segment Cox process. We 
see in this plot that "Shuffle" is able to recover the length of the 
longer segments more efficiently when samples are constructed 
out of 2D simulations instead of 3D simulations of segment Cox 
process. 
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Figure 11. This shows x 2 / u a t different shuffling lengths (L) 
for a mock SDSS NGP strip constructed out of 3D segment Cox 
process and 2D segment Cox process simulated with two different 
segment length I = 2fM _1 Mpc and I = SO/i^Mpc. The black 



dotted line indicates X l v = 1- 
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